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ABSTRACT 



We confirm the reality of a revers al of the sign of th e Faraday Rotation Mea- 
sure in the Galactic plane in Cygnus (ILazio et al Ill990l ). possibly associated with 
the Cygnus OBI association. The rotation measure changes by several hundred 
rad/m 2 over an angular scale of 2 — 5°. We show that a simple model of an 
expanding plasma shell with an enhanced density and magnetic field, consistent 
with observations of Ha emission in this part of sky, and physically associated 
with a superbubble of the Cygnus OBI association, can account for the magni- 
tude and angular scale of this feature. 



Subject headings: interstellar medium:bubbles — interstellar medium:HII regions- 
interstellar medium: magnetic fields 



1. Introduction 

Faraday rotation is one of the best diagnostics of the plasma properties of the inter 
stellar medium, in particular the st r ength and direction of the interstellar magnetic field 



(IClegg et al 1 11992 ; iFrick et al I l200ll ; lHaverkorn et al 1 12004 ; IVallee 1 12004 ; lHaverkorn et al 



20061 ). When radio emission from a source of polarized radio waves propagates through a 
plasma such as the interstellar medium, the measured polarization position angle \ is given 
by 



X = Xo + 



27rm 2 c 4 



nB ■ ds 



A 
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where xo is the intrinsic polarization position angle, assumed independent of wavelength (A), 
n is the plasma density, B is the vector magnetic field, and ds is an incremental pathlength 
interval along the line of sight. The quantities e, m e , and c are the usual fundamental 
physical constants of the elementary electrical charge, the mass of an electron, and the speed 
of light. The quantity within square brackets is termed the rotation measure (RM) and is the 
quantity which constitutes a plasma diagnostic of the interstellar medium (ISM). Equation 
(1) is in the cgs system of units. Conversion of the rotation measure to the more common 
SI units of radians/m 2 is done by multiplying the cgs value by a factor of 10 4 . 

The rotation measure is a signed quantity, depending on the density-weighted line of 
sight component of the magnetic field. The sign of the RM is of physical interest because 
it diagnoses the net polarity of the magnetic field along the line of sight. As may be seen 
in equation (1), an increase in \ with increasing A indicates that RM, and thus the path- 
averaged magnetic field, are positive. 

Measurements of rotation measure have been made along hundreds of lines of sight 
through the galaxy, using many radio telescopes. These observati ons have been used to de- 
termine the large scale structure of the Galactic mag netic field (e.g. Clegg et al 19921; 



1996 



Haverkorn et al 1120041 ) 



2004) as well as properties o f turbulence in the ISM jLazio et al Ill99pl : Minter and Spangler 



Vallee 



ar- 



This paper is particularly concerned with results reported by lLazio et al Nl990h . Pol 
ization measurements at three frequencies were made with the Very Large Array (VLA)0of 
8 extended radio sources observed through the Galactic plane in Cygnus, in the vicinity 
of the Cygnus OBI association. Previous radio observations had shown this region to be 
one of strong radio wave scatt ering, presumably due to p l asma turbulence associated with 



the Cygnus OBI associatio n (ISpangler and Cordes I Il988l ; iFey. Spangler. and Mutel I Il989 



Spangler and Cordes 19981). The distan ce to the Cygnus OBI association is estimated to be 
1.8 kpc iNichols-Bohlin and Fesen~lll993h . and thi s value for the dist ance will be employed in 
this paper. The main point of the observations of lLazio et al I (119901 ) was to determine prop- 
erties of the plasma turbulence on large spatial scales, manifest by differences ARM in the 
rotation measures on adjacent lines of sight to different parts of a background e xtragalactic 



radio source, and between lines of sight to different sources. The observations of lLazio et al 



(119901 ) were consistent with a Kolmogorov spectrum of irregularities in the plasma around 
the Cygnus OBI association, extending from spatial scales as small as 10 s — 10 9 cm, which 



1 The Very Large Array is an instrument of the National Radio Astronomy Observatory. The NRAO 
is a facility of the National Science Foundation, operated under cooperative agreement with Associated 
Universities, Inc. 
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are responsible for diffractive interstellar scattering, to as large as 10 17 — 10 18 cm, scales 
which cause Faraday rotation changes between adjacent lines of sight. 

The magnitude of Faraday Rotation in this part of the Galactic plane was found to 
be sufficiently large to cause difficulty in extrac ting RM va l ues fr om the observations at 
frequencies of 1.44, 1.65, and 4.88 GHz made by lLazio et al I (119901 ) . The basic problem is 
that for RMs of several hundred rad/m 2 , there can be several turns of the position angle 
between 1.65 and 4.88 GHz, and even a turn for more) between 1.44 and 1.65 GHz. This is 



illustrated graphically in Figure 3 of lLazio et al I (1l99Cj ). A technique was employed whereby 
the optimum set of values of RM consistent, in a least-squares sense, with the position angle 
measurements was retrieved. In the cases of four sources, 2004+369, 2005+368, 2007+365, 
and 2011+360, secure enough values for the average RM to the sour ce were obta i ned t o 
permit "unambiguous" rotation measures to be reported in Table 3 of lLazio et al I (jl990l ). 
The RM values reported were all positive, with values between 228 and 850 radians/m 2 . 



Clegg et al I (119921 ) presented polarization observations yielding RM for 56 extragalactic 



radio sources viewed through the Galactic plane, of which 4 were in the vicinity of the Cygnus 
OBI association. The RM values reported were large (several hu ndred radian s /m 2 ) , and 
negative. There w ere no sources in common between the samples of lLazio et al I (ll990l ) and 
Clegg et al I (Il992l ). This apparent large change in the RM (from several hundred radians/m 2 
positive to a similar magnitude, but negative) within a small angular distance suggeste d 
either (a) an error in the observations and/or data reduction procedure by lLazio et al I (ll990l ). 
so that the sign change did not exist, or (b) the existence of a remarkable plasma structure in 
the Galactic plane in this region, which could cause such a large, anomalous RM. We refer 
to this large change in the Galactic Faraday rotation measure over a small angular scale as a 
Rotation M easure Anomaly. A partial argument against possibility (a) was available even at 
the time of I Clegg et al (1992) . Processin g of the polarization position a ngle measur e ments 
made and published by lLazio et al I (119901 ) through the RM a lgorithm of IClegg et al I (119921 ) 
yielded the RM values given in Table 3 of lLazio et al I (Il990l ). 



The purp ose of this paper is to report observations made to confirm the large, positive 



RM values of lLazio et al I (Il990l ). Observations were made with the VLA at two, widely 
separated frequencies within the 5 GHz RF (radiofrequency) bandpass. Observations were 
ma de at this frequency because RM values of the magnitude reported by lLazio et al I (119901 ) 
and IClegg et al I (J1992J) would cause measurable rotation of several degrees of position angle 
or more between the two bands, but there would be no " n-7r ambiguity". As will be seen 



below (Section 3), these new observations do confirm the values reported in lLazio et al 



( I1990l ). indicating that a strongly Faraday- rotating structure is along the line of sight to 
these sources, and probably associated with the Cygnus OBI association. 
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The outline of this paper is as follows. In Section 2 we describe the polarimetric ob- 
servations made with the VLA at 5 GHz (C band), and optical observations made with the 
Wisconsin H-Alpha Mapper (WHAM) to investigate the plasma along this line of s i ght. I n 
Section 3, we compare the RM values obtained here with the values of lLazio et al I (119901 ). 
and see that they are consistent. In Section 4, we develop a simple model for the plasma 
structure of an expanding plasma shell, presumably a stellar bubble, associated with the 
winds and UV radiation of the Cygnus OBI association, and show that this model is capa- 
ble of reproducing the magnitude and angular scale of the observed RM anomaly. Finally, 
Section 5 summarizes and concludes. 



Observations 



2.1. VLA Polarimetric Observations 



Radio polarization measurements of four sources from the sample of lLazio et al I (jl990l ) 
were made on 13-14 December, 1999. The observations were made between 19:46 UT on 
the 13th to 3:42 UT on the 14 of December. The sources observed are given in Table 1. 
Column 1 g i yes th e J2000 name of the source, and column 2 gives the B1950 name used in 
Lazio et al I (119901 ). The remaining columns are described in Section 3 below. 



Three of the sources listed in Tab le 1, 2007+369, 200 9+367, and 2013+361, had rotation 
measures listed as "unambiguous" by lLazio et al I (jl990j). The fourth source, 2015+364, was 
not so listed but was observed in the present project because of the high polariz ed intensity. 



which facilitated a comparison of our observations and results with those of lLazio et al 
Jl990h . 



Observations were made at two frequencies, 4585 and 4885 MHz, within the C band RF 
bandpass. At both frequencies a 50 MHz intermediate frequency (IF) bandwidth was received 
and correlated. The 300 MHz separation between the two bands is about the maximum 
possible within the overall C band front-end bandpass, while still having similar gain and 
acceptable system temperature. For the purposes of this project, we consider polarization 
measurements at two frequencies as adequate for measuring Faraday rotation, rather than 
three frequencies as conventionally prescr ibed. As dis c ussed in the preceding section, given 
the rotation measure values published by lLazio et al I (Il990l ). the Faraday rotation between 
the closely-spaced frequencies of 4585 and 4885 MHz should be a small fraction of 180°. In 
this case, there is no possibility of the polarization position angle rotating through 180° or 
more between the two frequencies of observation (the u mr ambiguity"). At the time of the 
observations, the VLA was in the B array. 
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In addition to the four program sources given in Table 1, we also observed the source 
2007+404 as a phase and instrumental polarization calibrator, and 3C138 for purposes of 
calibrating the polarization position angles. The procedures of polarization calibration and 
mapping and analysis of the da ta are essent i ally i dentical to those in previous VLA polar- 
ization projects of ours, such as lLazio et al I (|1990l ). 



As in those previous investigations, the basic data products with which we worked were 
maps of the Stokes parameters Q and U, and maps of the polarization position angle x an d 
the polarized intensity L, 



X 



1 

- tan 
2 



-i 



(2) 



L=^(Q 2 + U 2 



Since the maps were made from 5 GHz data taken in the B array, the angular resolution of 
the maps was about 1.9 arcseconds (FWHM of synthesized beam). 

The procedure in measuring Faraday rotation from these observations was as follows. 
Maps of x an d L were made at both frequencies. The maps of x a t the two frequencies were 
subtracted from each other, giving a map of Ax, where the definition and significance of Ax 
are 



A X = x(4585) - x(4885) = RM [X 



4585 



^4885] 



(3) 



RM 



Ax 



X 2 

A 4585 



\ 2 

A 4885 



= 34.3Ax(o) rad/m 

In equation (3) x(4585) is the polarization position angle measured at a frequency of 4585 
MHz, A4585 is the corresponding wavelength at that frequency, and A%(o) is the difference 
in the polarization position angles at 4585 and 4885 MHz, in degrees. In the first two lines 
of equation (3), x an d Ax are in radians and A 458 5 and A4 885 are in meters, and RM has its 
conventional units of rad/m 2 . In the convenient formula of the third line of equation (3), 
RM is also in rad/m 2 , and Ax(°) is in degrees. 

Illustrations of our data reduction procedures are given in Figures 1 and 2, for the sources 
2007+369 and 2015+364, respectively. In each figure, panel (a) shows the polarization 



re radio sou r ce. A 
Lazio et al I f|l990h 
Jl990h 



Lazio et al 



position angle map, with variation reflecting the intrinsic structure of t 
comparison of these maps with the corresponding ones in Figure 2 of 
shows entirely satisfactory agreement, even though the observations of 
were made in the C array, and had roughly 1/3 the resolution of the present maps. This 
statement i s true for all four sources in common between the present project and that of 
Lazio et al I Jl99ph . 
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Fig. 1. — Polarization maps of the radio source 2007+369 at a frequency of 5 GHz. (a; 
left) Shows the map of the polarization position angle x at a frequency of 4.885 GHz. The 
contours are those of total intensity, and are given at -3,-1,1,3,5,10,20,30,50,70, and 90 % 
of the peak intensity, which is 3.54 mJy/beam. The restoring beam is shown in the lower 
left corner, and is 1.9 x 1.9 arcseconds. (b; right) Shows the distribution of position angle 
difference Ax between the position angles at 4585 and 4885 MHz 



The A% maps in panel (b) of both figures show a rotation between the two frequencies, 
which can be related to RM via equation (3). For both sources, Ax is nearly uniform 
across the source, being about 20° for 2007+369, and nearly zero for 2015+364 (see Table 
1 for exact values). The data in these figures clearly indicate that large Faraday rotation 
variations exist across the Cygnus OBI region. 

For each source we chose a number of regions (> 1) in which the polarized intensity 
was high enough for a reliable position angle measurement. In each region, we made a 
measurement of Ax at the position where L was locally a maximum. The errors were 
calculated on the basis of the conventional standard error in the polarization position angle 
measurement at a single frequency, a x = 

The results of these measurements are shown columns 3, 4, and 5 of Table 1. Column 3 
gives the source subcomponent for the measurement (e.g. Np for "north preceding", Sf for 
"south following", and cc for "central component", etc), column 4 gives the Ax measurement 
for that component, and column 5 gives the rotation measure to that component, calculated 
with equation (3). 
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Fig. 2. — Same as Figure 1, except for 2015+364. The peak intensity of the map is 7.11 
mJy/beam. 

2.2. WHAM Ha Observations 



We also analy sed data from the University of Wisconsin H-Alpha Mapper (WHAM; 
Haffner et al 1120031 ) for this part of the sky. The WHAM spectrometer samples the emission 
measure (EM = J n 2 e ds) of plasma along the line-of-sight at one-degree spatial resolution, 
allowing us to trace structures that could be responsible for the RM anomaly. We exam- 
ined both total intensity (Ha emission integrated over the whole line profile) as well as 20 
km/sec-wide channel maps centered from -60 km/sec to +60 km/sec with respect to the 
local standard of rest. 



3. Rotation Measures to Four Sources in Cygnus 

In this section we discuss the significance of the new measurements to the four sources 



i n Cyg nus, as well as the implications of the already-published measurements of lLazio et al 
(1l990[ ). 



3.1. Comparison with Observations from lLazio et al I (Il990l ) 



The RM valu es from this study (Column 5 of Table 1) are compared with those of 
Lazio et al I (119901 ) in Figure 3. All data from Table 1 of this paper are displayed, with the 
exception of the measurement of the central component of 2013+361, since the polarization 



- 8- 



of this component was not measured by lLazio et al I (ll990l ). The values plotted on the 



RM(Lazio) 



Fig. 3 . — Comparison of the RM values emergent from this paper with those of lLazio et al 



(1990) 



abscissa of Figure 3 are taken from Table 3 of lLazio et al I (jl990l ) except for 2015+364. Thi s 
was not listed as a source with an "unambiguous rotation measure" by lLazio et al I Jl990h. 
Given the small value of RM from our measurements, it was clear that the 



Lazio et al 



(Il990l ) measurements were not affected by t he "n n ambiguity", so we fit the position angle 
measurements from Table 2 of lLazio et al I (Il990l ) for the RM value. 



Figure 3 is the main observational result of this paper. The measurements from the 
current investigat ion are consistent in sign and magnitude with the values published by 
Lazio et al I (Il990l ). The apparent discrepancy for the Nf component of 2013+361 is a con- 
sequence of t he weak pola r ized i ntensity of this feature, and the correspondingly large error. 
Errors for the lLazio et al I (119901 ) measurements are plotted, but for most sources and source 
components are smaller than the plotted symbols. 

We confirm the large, positive rotation measures for the sources 2007+369 , 2009+367 



and 2 013+361. Given the accuracy of the RM values for these sources from lLazio et al 



( 119901 ). there is no reason to doubt the sign and magnitude of the fourth source with an 
unambiguous RM, 2004+369, which also had a large, positive RM (~ 820 rad/m 2 ). 

The new measurements also confirm the large spatial gradient in RM across this region, 
a change which can be seen to be at least 800 r ad/m 2 in ou r data alone, and approximately 
1300 rad/m 2 when the data sets of this study, lLazio et al I (119901 ) . and IClegg et al I (119921 ) 
are combined. 
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Finally, both our data and the refit position angle data from lLazio et al I (119901 ) show 

ablish a "handshake" between the results of 



Lazio et al 


(1990 


) and 


Cleee et al 


(1992) 



discuss its nature. 



4. The Source of the RM Anomaly 



In this section, we discuss the nature of the RM anomaly which has been confirmed 
by the present observations. We discuss a simple, physically based model for the plasma 
structure that could produce such a rotation measure anomaly, and also review the evidence 
for other such objects from observations in the literature. 

To begin the discussion, Figure 4 shows all RM data in the vicinity of the Cygnus OBI 
association. The rotation measure data are represented in the standard format of having 
solid circular symbols represent positive RM, open circles representing negative RM, and the 
plotted diameter of the circle being an indicator of the magnitude o f the R M. We show our 



measu rements, supplemented by the measurements of lLazio et al I (Il990l ) and IClegg et al 



(119921 ). In view of the large range of rotation measures in this part of the sky, the diameter 



of the plotted symbol in Figure 4 is proportional to the logarithm of the absolute magnitude 



Lazio et al 


(1990) and 


Cleeg et al 


(1992) 



The RM measurements are plotted on a background representing the intensity of the 
integrated Ha emission in this part of the sky from the WHAM instrument. The Ha emission 
shows a possible half arc of emission with an angular radius of 2.3°, centered on I = 75.0° 
and b = 1.0°. The lower half of this arc is obscured by heavy interstellar extinction in 
the "Aquila Rift". Radio continuum and recom b inatio n line observations of this region, 
summarized in Figure 6 of ISpangler and Cordes I (119981 ). show evidence for ionized gas at 
and below the Galactic equator at I = 76°, and extending to b = 1° at / = 74°. Interstellar 
extinction obscures the emission from this plasma in the WHAM image. In the discussion 
in Section 4.1 below, we will assume the physical reality of this shell, and attribute to it 
some properties of a stellar bubble. Although its existence and nature are by no means 
certain, this shell may correspond to the dense shell of shocked and photoionized interstellar 
gas associated with the bubble (or superbubble) produced by the Cygnus OBI association 
( IWeaver et al 1119771 ). A projection of the positions of the O and B stars in the Cygnus OBI 
association (taken from Humphreys 19781 ) on the Ha image shows that most of the stars fall 
within, or slightly beyond, the possible shell, which mildly corroborates our interpretation 
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of it as the swept-up plasma of an interstellar superbubble. This point was also made in 



Nichols-Bohlin and Fesen I (119931 ). 



4.1. Physical Model of the RM Anomaly 

In this section, we explore a simple physical model for the superbubble shell associated 
with Cygnus OBI, and see if it can produce a RM anomaly of the sort observed and displayed 
in Figure 4. We make no claims as to the uniqueness of this model. Rather, the goal is to 
show that structures of the sort known to exist in the vicinity of OB associations can produce 
a Faraday Rotation anomaly with a magnitude and angular extent of the sort we observe. 

Fig ure 5 shows a cartoon illustrating the bubble, which incorporates features from the 



model of I Weaver et al I (119771 ). The radius R defines the location of the outer shock. Inside 
Rq is a shell of thickness A, which contains dense, photoionized ISM material. It is this dense 
plasma which, in our model, produces the Ha shell seen in the WHAM image of Figure 4, 
and contributes to the enhanced rotation measure. 

The interstellar magnetic field (B ) will also be modified by this shell, and it contributes 
to the enhanced RM. We assume that it is uniform across the region exterior to the bubble, 
and is modified in the standard way by a magnetohydrodynamic (MHD) shock, in which the 
component normal to the shock front B n is unchanged, while the component in the shock 
plane B± is amplified by a factor X, Bj_ 2 — XB±i, where 1 and 2 refer to the upstream and 
downstream regions of the shock. In the case of a strong shock, X = 4. 

As indicated in Figure 5, the line of sight (LOS) perforates the shell at two points, I 
(ingress) and E (egress). The normal vectors to the shock at these two points are 

&i = %- (4) 

Kq 

where r} and Te are the vectors from the center of the shell to the ingress and egress points. 

The normal and perpendicular magnetic fields in the upstream regions at the ingress 
and egress points are 

Bni = {ni ■ S )nj (5) 
B nE = (n E ■ B )n E 

and 

Bxi = Bo~ Bni (6) 
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B±E — Bq — B n E 

The net magnetic fields in the shell (the downstream region) at the ingress and egress points 
are then given by 

Bn = B n i + XB ±I (7) 

— * — * — * 

Be2 — B n E + XB±e 

Finally, the z components of the magnetic fields in the shell at the ingress and egress points, 
B Iz and B Ez , which determine the Faraday rotation, are obtained by taking the dot product 
of equations (7) with the unit vector in the z direction (the direction from the source to the 
observer), e z . 

As seen from equation (1), the rotation measure is an integral along the line of sight. 
In a realistic model of a shell, the line of sight component of the magnetic field B z as well as 
the density would vary along the line of sight through the shell. Given the simplified nature 
of the model presented here, we approximate the RM as 

RM(0 = [B zI + B zE \ (8) 

where C is the collection of atomic constants within the curved brackets in equation (1), 
which has the value 2.631 x 10~ 17 in cgs units, or 0.81 if L is in parsecs and B z is in 
microGauss. The variable n is the assumed uniform plasma density in the shell, B z j and 
B zE are the downstream (postshock), line of sight components of the magnetic field at the 
ingress and egress points, respectively. The length of the chord through the shell is denoted 
by where £ is the distance of the line of sight from the center of the shell, as illustrated 
in Figure 5. Equation (8) assumes that all plasma characteristics are uniform in two branches 
of the line of sight, each of thickness L/2 as shown in Figure 5. The first branch is from 
the ingress point to an interior point, and the second branch extends from the interior point 
to the egress point. The values of these interior points depend on whether the parameter 
£ < R 1 = Rq — A, or £ > Ri. In the former case, the interior points are points where the 
line of sight intersects the inner surface of the shell. In the latter case, both interior points 
are the same as the midpoint of the chord through the shell. This definition is clarified by 
reference to Figure 5. Equation (8) is the fundamental equation in our analysis. 

It may easily be shown that the rotation measure through a shell of the sort shown in 
Figure 5 is of order RM = 2CnB z0 R . The maximum value of the rotation measure occurs 
for chords with £ > R-y, and is ^RM for a compression factor X = 4. The difference ARM 
between the maximum rotation measure and the RM of the central chord depends on the 
ratio but is of order ARM ~ |i?M . 
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4-1.1. Estimating the Plasma Density in the Shell 

Equation (8) takes as its input several properties of the plasma shell illustrated in the 
cartoon of Figure 5, and perhaps seen in the Ha map of Figure 4. Obviously, one of the 
most important is the assumed constant plasma density n in the shell. This parameter can 
be estimated from the intensity of the Ha emission shown in Figure 4. The re lation between 



the m easured Ha intensity in Rayleighs, I Ha and the emission measure is (IMadsen et al 



2006|) 



EM = 2.75T 4 a % a (9) 

where EM is in the customary units of cm _6 -pc, and T 4 is the plasma temperature in units 
of 10 4 K. If we adopt a customary value of T 4 = 0.8, this relation becomes 

EM = 2.2hI Ha (10) 

We measured the In a values from the WHAM map for each of the lines of sight through the 
shell to a radio source. The resultant emission measures showed modest variation from one 
line of sight to another, with a total range of 192 - 360 cm~ 6 -pc. 

It is obvious that interstellar extinction will reduce these emission measures below their 
true values, with corresponding underestimates of the plasma density in the shell. We 
estimated the roll of extinction in the following way. 



Spangler and Cordes I (119981 ) analysed a number of lines of sight through the same part 



of the Cygnus OBI ass ociation (compare the locat ions of the sources in Figure 4 of this 



paper with Figure 6 of ISpangler and Cordes Ill998l ). in that case studying the magnitude 



of interstellar scattering. As part of that analysis, values for the emission measure along a 
number of lines of sight to extragalactic radio sources were determined from radio continuum 
brightness temperatures, which are not subject to interstellar extincti on. The values for the 



emissi on measure along those lines of sight are given in Table 2 of ISpangler and Cordes 



( 119981 ) . We measured the Ha intensity along these same lines of sight, and used them to 



obtain an estimate of the emission measure from the WHAM data, whi ch was compared with 



the ra dio-based emission measures. For the five sources in Table 2 of ISpangler and Cordes 



( 119981 ). the ratio a of the radio based EM to the optically-based one ranged from 1.4 to 4.2, 
with a mean of 2.6. The range in values of a may well represent spatial variations in the ex- 
tinction from one line of sight to another. Furthermore, since the beam of the radio telescope 
used for the radio continuum measurements (Max Planck Institut fur Radioastronomie at 
Effelsburg) and WHAM are different, the effective extinction in the two beam solid angles 
may not be the same. Nonetheless, to obtain a first-order estimate of the effects of extinction 
in this region, we corrected the WHAM emission measures by the mean value of a = 2.6. 
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Given all of this, our estimate for the density of the shell, calculated along each line of 
sight, is given by 

l2.25aI Ha 

The density calculated in equation (11) assumes the shell is a spherical annulus with 
uniform density. If the plasma is clumped with a filling factor less than unity, the density 
calculated in (11) will obviously be less than the density in the clumps. The rotation measure 
calculated under the assumption of a uniform shell with density given by (11) will exceed the 
true rotation measure from a shell with a clumped density. It should be noted from equation 
(11) that the derived density is weakly dependent (proportional to the square root) on the 
imperfectly-known parameter a. 

For each line of sight, the source-specific value of £ was used. Calculation of also 
requires estimates of the distance to Cygnus OBI, which we take to be 1.8 kpc, and the 
thickness of the shell. From the WHAM data shown in Figure 4, we estimated Rq = 71 pc 
and A = 39 pc. The values for the density which resulted spanned a fairly narrow range of 
2.3-3.9 cm" 3 . 

Estimates of the external magnetic field Bq were made as follows. We assumed that the 
magnitude of the ISM field i s |j? n l = 3 — 5/iG, as indicated by numerous independent studies 



using Faraday rotation (e.g. iMinter and Spangler Ill996t lHaverkorn et al 1 12004 ; lHaverkorn 



20071 ). We adopted a value of B = 4.0/xG. The orientation of the interstellar magnetic field 
with respect to the line of sight at the location of the Cygnus OBI association bubble is 
unknown, and was left as an adjustable free parameter. 

To model the rotation measures for the sample of sources shown in Figure 4, we adopted 
values of Rq, A, n, and B as described above. Additional free parameters which affect 
the RMs are the angle 9 between the interstellar magnetic field and the line of sight at 
the location of the Cygnus OBI association, and an offset constant RM Q ff, which is the 
background rotation measure in this part of the Galactic plane. Values for the latter two 
parameters were chosen which gave a good match to the observed rotation measures. The 
calculated model RMs used the specific value of n for each line of sight as obtained from 
equation (11). The model RMs thus incorporated available observational information, rather 
than use a single mean value for n as would be required by strict adherence to the model. 

The parameters of this model for the shell associated with the Cygnus OBI superbubble 
are given in Table 2 as "Model A". The rotation measures emergent from this model are 
shown in Figure 6, in the same format as Figure 4. Figure 6 is qualitatively similar in appear- 
ance to the true measurements, given in Figure 4. The range of model RMs is from +665 to 
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359 rad/m 2 . We als o carri ed out calculations for a model shell similar to that discussed by 



Spangler and Cordes I (119981 ). and based on observations from the literature and theoretical 
arguments. This second model ("Model B") has Rq = 94 pc and A = 13 pc. Model B has 
it own set of associated shell densities obtained from the emission measure measurements, 
and the free parameters RM Q ff and 9 were adjusted to give a good representation of the 
observations. This model also satisfactorily represented the observed range of RMs, with 
the model RM values ranging from +784 to -274 rad/m 2 . The parameters for Model B are 
also given in Table 2. 

The ability of these simplified shell models (with reasonably accurate values of the shell 
parameters) to reproduce the magnitude and angular scale of the RM values we observe is 
an important scientific result of this paper. We do not claim that the models presented in 
Table 2 are unique. It is obvious from equation (8) that the plasma density, the magnitude 
of the upstream magnetic field, and the angle between the interstellar field and the line of 
sight are highly correlated. 

The calculations presented above are sufficient to demonstrate that a physical structure 
similar to the plasma shell associated with the Cygnus OBI superbubble would produce a 
rotation measure anomaly similar to that which we observe. A more realistic model, incorpo- 
rating additi onal independent information on the properties of the Cygnus OBI association 
superbubble ( jSaken et al lll992l ; lNichols-Bohlin and Fesen Hl993uSpangler and Cordes 1l99Sh 
would probably result in better agreement between model observables and observations. 



4.2. Comparison with Other Anomalies Reported in the Literature 



In this subsection, we discuss the evidence for other, similar features in the sky, which 
might be physically related to the plasma shell around Cygnus OBI. Jacques Va ll ee ha s 



for many years reported and discussed RM features on the sky (e.g. iVallee I Il993l . 120041 ). 
which he refers to as interstellar magnetic bubbles. In comparison with the object discussed 
here, the interstellar magnetic bubbles are larger in angular extent, ar e at higher Galactic 
latitudes, and induce smaller changes in the background RM. Alt hough I Vallee I (119931 . 120041 ) 
does not model these variations in the same way as we do, he does posit a connection with 
stellar OB associations. Further consideration of the similarity (and difference) of the objects 
discussed by Vallee, and the plasma shell associated with the Cygnus OBI association would 
be worthwhile. 

An object which is in some ways more similar to that described in this pape r is th e 
rotation measure anomaly in the Galactic plane at I ~ 92°, discovered by lClegg et al I (119921 ). 
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Table 1. Sources Observed and Faraday Rotation Measurements 



Source (J2000) 


Source (B1950) 


Comp. 


Ax(o) 


RM (rad/m 2 ) 


2007+369 


2005+368 


Np 


20.1 + 0.7 


689 ± 24 






Sf 


21.4 + 0.4 


732 ± 12 


2009+367 


2007+365 


Sf 


7.8 + 1.1 


269 ± 40 






Np 


5.9 + 1.4 


201 ± 48 


2013+361 


2011+360 


Sp 


4.5 + 0.9 


154 ± 30 






cc 


10.9 + 6.2 


373 ± 213 






Nf 


-1.1 + 6.6 


-37 ± 226 


2015+364 


2013+362 


Sp 


-0.2 + 0.2 


-7.5 + 6.2 






Nf 


-3.3 + 0.4 


-114 + 15 



Table 2. Models for Shell Associated with Cygnus OBI association 



Shell Parameter 


Model A 


Model B 


Ro (parsecs) 


71 


94 


A (parsecs) 


39 


13 


n (cm -3 ) 


2.3-3.9 


4.3-5.3 


B ( M G) 


4.0 


4.0 


RM off (rad/m 2 ) 


-723 


-676 


6 (degrees) 


59 
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In that paper IClegg et al I (119921 ) suggest that it may be similar to our anomaly in Cygnus. 
Th e I g 92° anomaly has bee n mapped in more detail, and with greater density of sources, 
by iBrown and Taylor I (120011 ). who identified it as one of three similar anomalous regions 
(not including Cygnus OBI), although the feature at / ~ 92° is the clearest case in their 
sample. The change in RM due to this region is of the order of several hundred to 1000 
rad/m 2 , which is quite comparable to the Cygnus OBI feature. The angular extent is also 
similar, with an angular diameter of 2.5° (IBrown and Taylor 1 120011 ). However, in spite of 
these similarities between the two features, there are no obvious HII regions or interstellar 
bubbles similar to Cygnus OBI along the line of sight to the I ~ 92° anomaly. Although 
further attention is clearly necessary, this simple, qualitative astronomical observation would 
seem to preclude a plasma structure or stellar bubble of the sort we have discussed in Section 
4.1. It may be the case that different t ypes of object s can produce such RM anomalies, a 
conjecture supported by the catalog of I Vallee I (120041 ) . Given the similarity in the Faraday 
rotation of the Cygnus OBI and I ~ 92° features, we cannot exclude the possibility that the 
bubble or plasma shell associated with the Cygnus OBI association is not responsible for 
the strong Faraday rotation there. 

The best approach for further illuminating this issue would be to examine OB associ- 
ations in the direction of the Galactic anticenter, where confusion of multiple objects along 
the line of sight is not an issue. Such an investigation could include both quantitative mea- 
surement of the properties of the Ha-emitting shells, as well as measurements of, or upper 
limits to, Faraday rotation anomalies. 



5. Summary and Conclusions 



The conclusions of this paper are as follows. 



1. Polarization observations at the closely spaced frequencies of 4585 and 4885 MHz with 
the VLA have confirmed the large, positive Farad ay rotation measure s for the sources 
2007+369, 2009+367, and 2013+361 reported by lLazio et al I (|l99oh on the basis of 
observations at very different frequencies, and analysed in different ways. 



2. The observations in this paper, together with those of lLazio et al I (Il990l ) and lClegg et al 
(119921 ). appear to show the presence of a "Faraday Rotation Anomaly", in which the 
RM changes by several hundred radians/m 2 , or more, over an angular scale of 2-4 
degrees. The RM changes by about 800 radian/m 2 between the sources observed in 
this proje ct, separated by an ang ular distance of ~ 2°. Inclusion of other sources ob- 



served by lLazio et al I (119901 ) and IClegg et al I (119921 ) leads to differences as large as 
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1300 radian/m 2 over angular separations which are only slightly larger. 

3. This Rotation Measure Anomaly occurs for sources which are viewed through a part 
of the Galactic plane which contains the Cygnus OBI association. Ha images from 
the WHAM telescope show what may be a plasma shell caused by the stellar bubble 
or superbubble produced by the association. 

4. We have developed a simple physical model of a plasma shell associated with a stellar 
bubble which appears to be present in Cygnus OBI. Using independent estimates of 
the shell size, thickness, and mean density, as well as the strength of the interstellar 
magnetic field outside the shell, we have produced RM gradients which are consistent 
with the observations. 

5. The physical relation of this object to other, somewhat similar or very similar objects 
reported in the literature is not entirely clear, and would benefit from further research. 

This research was supported at the University of Iowa by the National Science Foundation 
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Fig. 4. — Faraday Rotation in the vicinity of the Cygnus OBI association. Positive and 
negative RM are indicated by solid and open circles, respectively. The diameter is propor- 
tional to the logarithm of the absolute ma gnitude of the RM . The data are fro m the present 
paper (4 sources), the source B2004+369 (ILazio et al Ill990l ). and 4 sources from lClegg et al 



(11992). Black symbols p ort ray measurements from this paper, and white symbols are data 



from ILazio et al I (119901 ) or Oegg et al I (119921 ) . The shaded background illustrates the in- 
tensity of the Ha emission as mapped by WHAM. White corresponds to higher intensity, 
as indicated by the grayscale bar at the right. The peak intensity is 197 Rayleighs. The 
WHAM data show an approximate arc of emission with a radius of 2.3° centered on Galac- 
tic coordinates I = 75.0° and b = 1.0°. The lower half of the possible arc is absorbed by 
interstellar extinction in the "Aquila Rift" . 



Fig. 5. — Cartoon of a simplified model of a plasma shell associated with a stellar bub- 
ble produced by the Cygnus OBI association. The cartoon serves to define several of the 
variables which are used in the discussion in Section 4.1. 
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Fig. 6. — A set of model rotation measures along the observed lines of sight, in the same 
format as Figure 4. For this calculation, we used plasma shell densities of 2.3-3.9 cm -3 , and 
an external interstellar magnetic field of 4.0/iG oriented at an angle of 59° with respect to 
the line of sight. 



